SIGNATURES OF LONG-LIVED SPIRAL PATTERNS. 
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ABSTRACT 

Azimuthal age/color gradients across spiral arms are a signature of long-lived spirals. From 
a sample of 19 normal (or weakly barred) spirals where we have previously found azimuthal 
age/color gradient candidates, 13 objects were further selected if a two-armed grand-design pat- 
tern survived in a surface density stellar mass map. Mass maps were obtained from optical and 
near-infrared imaging, by comparing with a Monte Carlo library of stellar population synthesis 
models that allowed us to obtain the mass-to- light ratio in the J-band, (M/L)j, as a function of 
(g — i) versus (i — J) color. The selected spirals were analyzed with Fourier methods in search 
for other signatures of long-lived modes related to the gradients, such as the gradient divergence 
toward corotation, and the behavior of the phase angle of the two-armed spiral in different wave- 
bands, as expected from theory. The results show additional signatures of long-lived spirals in at 
least 50% of the objects. 



Subject headings: galaxies: kinematics and dynamics 
galaxies: stellar content — galaxies: structure 



galaxies: photometry — galaxies: spiral 
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1. Introduction. 

Understanding the origin of spiral arms, bars, 
and rings is key for discerning the long-term evolu- 
tion of disk galaxies. The density wave (DW) the- 



ory, originally conceived by Lindblad (1963) and 
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developed in a linear form by Lin & Shu (19641, 



has been commonly brought up to explain spiral 
structure. One fundamental prediction of this the- 
ory is the presence of offsets between tracers of 
the different stages of star formation and evolu- 
tion ( Roberts|1969 ). A manifestation of these off- 
sets would be age/color gradients across arms, as 
a result of the interaction between a spiral pattern 
with an angular speed that is nearly constant with 
radius, and stars and gas with differential rotation. 
The gas overtakes the spiral pattern inside coro- 
tation (CR, i.e., the radius where the pattern and 
orbital angular velocities are the same), and is left 
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behind outside CR. As gas feels the gravitational 



potential near the spiral arms, it shocks (Roberts 



1969 Gittins & Clarke 2004). This produces the 



piling up of dust and molecular gas that is com- 
monly found in the concave side of the arms, if the 
process occurs inside CR and spiral arms trail. If 
the shock triggers star formation, an age (or color) 
gradient is expected across spiral arms as the new 
stars age and drift away from the arms. Young 
stars would be located immediately downstream 
the dust lanes, followed by increasingly older stars 
in the direction of rotation. In part due to the dif- 
ficulties encountered to detect such gradients, the 
relevance of spiral DWs for disk galaxies has been 



ct al. 



ct al. 



(2012 



repeatedly challenged. Recent examples are Foyle 
20111 ), |Grand et al.| p012a|b| ), and |Ferreras 



In a series of papers, we have detected and 
studied azimuthal color/age gradients across spi- 
ral arms and bars, thus establishing a relation 
between spiral or bar dynamics and star forma- 



tion ( 


Gonzalez & Graham 


1996 Martmez-Garcfa 


et al. 2009a|b 


Martfnez-Garcfa & Gonzalez- 


Lopezlira 20111. 


Our studies suggest that az- 



imuthal gradients can be found in almost every 
galaxy, but are difficult to detect and do not ap- 
pear widespread within a single galaxy. |Efremov 
(2010) has noticed an inverse correlation between 



the signatures of a spiral shock, and the presence 
of chains of star complexes ("beads on a string") 
in spiral arms; he has also observed that spiral 
shocks seem to be associated with an irregular 
magnetic field (Dobbs & Price 2008), whereas the 
"beads on a string" configuration correlates with 
the existence of a regular magnetic fieldj^j 

The "infant mortality" of star clusters ( |Lada k. 
Lada||2003| |Goodwin fc Bastian||2006[ ) , that leads 
to their disruption in the first 10-30 Myr, may pre- 
vent the formation of age gradients even if shocks 
trigger star formation. Stellar winds and super- 
novae may drive out the gas not used in star for- 
mation, causing a change in the gravitational po- 
tential. This would lead to the disintegration of 
the clusters, as stars become unbound. Coher- 
ent age gradients must first survive the "infant 



1 On smaller scales than color gradients, Herschel Space Ob- 
servatory | |Pilbratt et al.|2010[ ) has shown that stars form 
along gaseous filaments, whose width is surprisingly con- 
stant, ~ 0.1 pc, irrespective of central column density \Ar-\ 

|zoumanian et al.|201l[ ). 



mortality" process; then, after i age > 50 Myr, the 
"dissolution of stellar groups" scenario proposed 



by Wielen (1977) can also take place near spiral 
arms. These last process can explain the appar- 
ent "downstream decline" of the gradients already 



noticed by us (cf. Martmez-Garcfa et al. 2009a 



Martmez-Garcfa & Gonzalez-Lopczlira 20111. 
1.1. Pattern speed variation with Radius 



Recently, based on the analysis with the radial 
Tremaine- Weinberg (RTW) method of CO and HI 
data of several galaxies, it has been proposed that 
the spiral pattern speed fl p may actually increase 



with decreasing radius in some objects (Merri- 


ficld et al. 2006| |Meidt et al. 


2009| |Speights & 


Westpfahl|2012 1. According to 


Wada et al. (2011) 



and Grand et al. ( 2012a|b l, this behavior is also 
seen in simulations^ If this is indeed the case, 
there would be no real long-lived pattern with solid 
body rotation, and hence no significant offsets be- 
tween the tracers of the different stages of star 
formation would be expected (see, for example, 
Grand et al. 2012b). In particular, azimuthal color 



gradients across the arms would not be observed. 
In Martmez-Garcfa et al. (2009b I, we investi- 



gated the effect of having (and neglecting) non- 
circular streaming motions, that is, of using a 
purely circular dynamic model to derive spiral pat- 
tern speeds from observations of color gradients. 
Summarizing, based on semi-analytical solutions 
and MHD simulations, gas orbits were obtained for 
models with a fixed pattern speed. Stellar popula- 
tion synthesis (SPS) models of age gradients were 
incorporated to these gas orbits, and synthetic 
photometric observations were performed]^] The 
method of Gonzalez & Graham (1996 hereafter 



GG96) was then applied to measure f2 p . In spite of 
having a real pattern speed, constant with radius, 
higher pattern speeds were measured at smaller 
radii (see Fig ure [T|) . As mentioned in Martinez 



Garcia et al. (2009b), the reason for this bias is 
that gas streamlines in a steady rotating spiral 
shock turn somewhat along the arm after pass- 
ing through the shock. Hence, stars take longer 



Roca-Fabrega et al. (2013) find that a radially varying Q p 
occurs only in simulations with flocculent, transient arms. 
3 We assumed that young stars retain the velocity compo- 
nents of their parental molecular clouds. 
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to move away from the arm, and the observer will 
think there is a smaller difference between the pat- 
tern speed and the orbital frequency. f2 p will seem 
to follow f2, just like in the previously mentioned 
applications of the radial RTW method, except 
that in our case we can be sure that the trend 
is caused by the streaming motions of the young 
stars, rather than by a radial variation of f2 p . The 
bias is also found in real observations, as shown 
in 



ibrated with the same standard procedures used 



Martmez-Garci'a et al. (2009b) 



In an ideal scenario where f2 p is indeed constant 
for all radii and spiral arms trigger star formation 
in all their extension, color gradients are predicted 
to be located all along the arms (except near CR) . 
In reality, age gradients are detected only in a few 
small regions per galaxy, but they do seem linked 



to disk dynamics (GG96; Martmez-Garcia et al 



2009a b ; Martme z-Garcia &: Gonza lcz-Lopczlira 
2011 ). So far, we have detected gradients through 



visual inspection. In order to locate gradients in 
a more objective way, and to attempt to identify 
some that do not meet the theoretical expecta- 
tions of coherence or smoothness, in what follows 
we will adopt Fourier techniques to study normal 
(or weakly barred) spirals. 

2. Photometric data. 

Our data consist of deep photometric images of 
the 13 SA and SAB spiral^jdescribed in Martinez 
Garcia et al. ( 2009a[ ) . The data include mosaics in 
the optical g, r, and i bands dThuan & Gunn|1976 



Wade et al. 1979 



sec Table nT the near-infrared 
(NIR) J band, and the NIR H, K, K s or K' fil- 
ters. We have now added six new objects Q whose 
observation log is shown in Table [2j NGC 1703 
and NGC 6907, classified as SB by RC3, were in- 
cluded in the analysis because (apparently) they 
have no prominent bar structure, as appreciated 
in the images. The data were reduced and cal- 



4 NGC 578, NGC 918, NGC 1417, NGC 1421, NGC 3162, 
NGC 3938, NGC 4254, NGC 4939, NGC 5371, NGC 6951, 
N GC 7125, NGC 7126, and N GC 7753. As argued 
in |Ma rtmez-Gar cfa et al,| (2009a) , NGC 578 is peculiar, 
in the sense that it is the only object in the sample whose 
arms may end at CR and not at the outer Lindblad reso- 
nance (OLR). We also caution the reader that NGC 1421 
is seen nearly edge-on, with an inclination angle a ~ 76°. 

5 NGC 1703 (type SB), NGC 3001 (SAB), NGC 3338 (SA), 
NGC 4603 (SA), NGC 6907 (SB), and NGC 7083 (SA). All 
galaxy types are taken fro m the Third Reference Cata logue 
of Bright Galaxies (RC3, de Vaucouleurs et al.|1991 i. 



in Martmez-Garcia et al. 


( 2009a 


1, and 


Martinez- 


Garcia & Gonzalez-Lopezlira 


(2011 


). The optical 



(g, r, and i) calibration was done in the Thuan- 



Gunn system (Thuan & Gunn 1976 Wade et al 



1979 ) . The zero point of this photometric system 
is chosen such that the standard star BD+17°4708 
has g = r = i = 9.5 mag. The NIR J data were 
calibrated with images from the Two Micron All 
Sky Survey (2MASS; [Skrutskie et al.|[2006 |. 



3. Structural type of the spiral arms 

The DW theory predicts the existence of age 
gradients across spiral arms. Multi-armed or floc- 
culent galaxies, however, are not described by this 
theory. The short arms of the latter are probably 
unbound star complexes sheared by galactic differ- 
ential rotation, whereas the "knee-like'j^] shaped 
arms of MIOl-type galaxies might be formed by 
the gravitational instability of galactic disks; if 



so, they are transient features (Clarke & Gittins 
2006||Dobbs fc Bonnell|2008[ ) . The azimuthal pro- 
files (intensity versus 9) of such arms should be 
symmetric about the gravitational potential well, 
where the gas is yet to be transformed into stars; 
this configuration is indeed sometimes observed in 
certain galaxies (for example, in NGC 4535 and 
M61, a narrow dark lane is seen right at the cen- 
ter of a stellar arm). 

To analyze the validity of the DW theory, one 
should study galaxies whose structure is likely to 
be explained by it. Mainly, from the point of view 
of their spiral arms, there are three different types 
of galaxies (e.g., Efremov 2011[ ). These are (1) 
symmetric grand-design spirals (e.g., M81), with 
mass arms due to DWs; (2) multi-armed or "knee- 
like" spirals (e.g., M101); and (3) flocculent spiral 
galaxies (e.g., NGC 2841). However, even in the 
same galaxy, different arms may be formed by dif- 
ferent mechanisms. Thus, in the classical multi- 
armed galaxy M101 the inner arms are rather reg- 
ular, and in our own Galaxy the quite symmetri- 
cal inner arms transform into a polygonal multi- 



armed structure in the outer regions (Efremov 



2011). Albeit flocculent and multi-armed struc- 
tures have a low mass-to- light ratio (M/L) and are 
seen mainly in the optical (see, for example, the 



6 Consisting of straight fragments that join at angles of 
roughly 135°. 
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Table 1 
Filter Characteristics 



Filter 


A e ff 


FWHM 


S 


5000AA 


830AA 


r 


6800AA 


1330AA 


i 


7800AA 


1420AA 


J 


1.25/im 


0.29/jm 


Note 


. — Columns 


2 and 3: el- 


fcctivc 


wavelengths 


and widths, 



respectively, of the filters mainly 
used for analysis. 



Table 2 

Observation Log and Galaxy Parameters 



Object 


Filter 


Exposure (s) 


Telescope 


Date (month/year) 


* Parameters 


NGC 1703 


9 


5100 


CTIO 0.9 m 


3/94, 3/95 


Type: SB(r)b 




r 


4200 






PA (deg): 104 a 




i 


5400 






a (dog): 27.0 ± 12.1 




J 


2717 


CTIO 1.5 m 


2/94 


v zot (km s- 1 ): 55 ± 20 




K 


1281 






Dist (Mpc): 20.4 ± 1.7 




H 


2758 








NGC 3001 


9 


4500 


CTIO 0.9 m 


3/94, 3/95 


Type: SAB(rs)bc 




r 


5100 






PA (deg): 6 




i 


4500 






a (deg): 47.5 ± 4.9 




J 


2718 


CTIO 1.5 m 


2/94 


u rot (km s _1 ): 248 ± 20 




H 


2700 






Dist (Mpc): 35.6 ± 3.0 


NGC 3338 


a 


3600 


Lick 1 m 


4/94, 11/94 


Type: SA(s)c 




r 


4500 






PA (deg): 100 




i 


4200 




2/94, 4/94, 11/94 


a (deg): 51.9 ± 2.1 




J 


1200 


Lick 1 m 


12/94, 2/95 


w rot (km s" 1 ): 194 ± 9 




J 


1002 


Kitt Peak 1.3 m 


3/94, 11/94 


Dist (Mpc): 23.7 ± 2.1 




K s 


926 








NGC 4603 


a 


5400 


CTIO 0.9 m 


3/94, 3/95 


Type: SA(s)c 




r 


4500 






PA (deg): 27 




i 


4200 






a (deg): 43.6 ± 5.7 




J 


1393 


CTIO 1.5 m 


2/94 


v rot (km s" 1 ): 242 ± 32 




H 


1399 






Dist (Mpc): 29.4 ± 2.6 


NGC 6907 


a 


3600 


CTIO 0.9 m 


9/94 


Type: SB(s)bc 




r 


3900 






PA (deg): 46 




i 


3600 






a (dog): 35.6 ± 3.7 




J 


570 


CTIO 1.5 m 


9/94, 9/95 


■u rot (km s" 1 ): 244 ± 41 




K a 


355 






Dist (Mpc): 47.7 ± 4.1 


NGC 7083 


a 


3600 


CTIO 0.9 m 


9/94 


Type: SA(s)bc 




r 


3900 






PA (deg): 5 






3900 






a (deg): 52.9 ± 3.0 




J 


840 


CTIO 1.5 m 


9/94 


v [ot (km s _1 ): 226 ± 10 






392 






Dist (Mpc): 45.4 ± 3.9 



Note. — Column 6-galaxy parameters: Hubble types from the RC3; P. A.: position angles from RC3; inclination 
angle, a — cos -1 (b/a), where a/b is the isophotal diamet er ratio derived from the i?25 parameter in RC3; v Tat : 
galactic rotation velocity obtained from the HI data of |Paturel et al.| < |20(33| , corrected for inclination; dist : 
Hubble distance obtained from the RC3 heliocentric radial velocity and the infall model of |Mould et a!~| < |2000^ 
H = 71 ± 6 km s _1 Mpc -1 . 

a |Paturel et al.|f2000( . 
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classical example of NGC 309 in Block et al. 1994), 
in order to determine whether the arms in the NIR 
are indeed mass DWs it is important to disentan- 
gle the contribution of young stars and clusters at 
these longer wavelengths; such contribution can 
reach up to 20%-30% (e.g., Rix & Rieke 1993; 
GG96; Rhoads 1998; James & Seigar 1999; Patsis 
et al. 2001; Grosb0l et al. 2006). With the aim of 
determining the actual structural type of the spi- 
ral arms, resolved stellar mass maps are required. 
Mass maps get around the problem that the M/L, 
at any wavelength, is not constant across the disk 
and near spiral arms. 

In order to build resolved maps of stellar mass 
of our sample galaxies, we use the method of |Zi-| 



betti et al. (2009). The method relies on a Monte 



Carlo library of 50,000 stellar population spectra, 
constructed from the SPS models of IBruzual~fcl 



Chariot (2003), and Chariot & Bruzual (2007, pri- 



vate communication) . Each spectrum is computed 
by randomly drawing the model parameters (star 
formation history, metallicity, and dust attenua- 
tion^ ]) from adequate physical distributions (see 



also |da Cunha, et al.||2008[). We use the "MAG- 
PHYS" codf | of |da Cunha et al.| ( |2008[ ) to ex- 
tract colors from the spectral library. Each ele- 
ment in a theoretical two-dimensional color-color 
plot can be produced by several combinations of 
model parameters. The median M/L of all these 
combinations is then found, and attributed to the 
data with the same position in the observed two- 
dimensional color-color plot. Resolution elements 
in two-dimensional color space (in our case, g — i 
versus i — J) have a size of 0.05 x 0.05 mag 2 , in 
order to account for the observational error. We 
compute the median M/L in the J-band (see Fig- 
ure 0), and compare the resulting distributions 
with the colors obtained from the observations^] 



7 Treated as prescribed in the two-component dust model 
of |Charlot k. Fall| ( [2>jo7j] ) . 

8 http: / /www. iap.fr/magphys/magphys/MAGPHYS. html. 

9 All images are treated with the adaptive smooth- 
ing code "A DAPTSMOOT H" , developed by Ste- 
fano Zibetti (Zibctti 2009), and available at the 
URL http:/ /www. arcetri. astro, it /~zibetti/ Software / 

ADAPTSMOOTH.html The code enhances the S/N with 
a minimum loss of effective resolution, and keeps the 
photometric fluxes unaltered. The models are reddened 
before comparing with the observations, to take into 
account foreground Galactic extinction (Schlcgcl et al. 
|1998| . 



on a pixel-by-pixel basis. This comparison pro- 
duces a resolved map of (M/L)j. Finally, a rela- 
tive "mass map" is obtained, equal to the product 
of the (M/L) j map by the J-band image. Ab- 
solute mass maps can be worked out, taking into 
account the distance to the galaxies, but are not 
necessary for our purposes. 

The results are shown in Figures [3] and |4j where 
we compare the g and J frames with the mass 
maps at the same relative scale. We then proceed 
to select the objects in our sample that retain a 
spiral structure that may be explained by the DW 
theory. From this wsitap*] inspection we reject the 
following objects prior to the Fourier analysis: 

1. NGC 3162. There is no two-armed spiral 
structure appreciable in the mass map. 

2. NGC 3938. The mass map has a very low 
quality, due to a low signal-to-noise ratio 
(S/N) of the g-band mosaic. 

3. NGC 5371. Ring structures rather than spi- 
rals are discerned in the mass map. 

4. NGC 7083. The spiral structure looks floc- 
culent in the mass map. 

5. NGC 7126. No clear spirals can be appreci- 
ated in the mass map. 

6. NGC 7753. Spiral structure is not observed 
in the mass map. 

The remaining 13 objects were analyzed as de- 
scribed below. 

4. Fourier analysis of the gradients. 



4.1. Azimuthal Q(rJgi) Index Analysis. 



The hrst step in our Fourier analysis consists in 
generating Q(rJgi) mosaics, following the method 
of GG96. Q is defined as follows: 



Q{rJgi) = {r-J)- E }j J) A g-i), (1) 
E(g-i) 



10 Future works will include a more objective selection pro- 
cedure. 
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where E Sj~ J ) is the color excess term. This pho- 

E(g-i) _ _ r 

tometric index is reddening-free for screen absorp- 
tion, and also for a mixture of dust and stars, as 
long as Ty < 2. At the same time, Q(rJgi) is sen- 
sitive to supergiant stars, and therefore traces star 
formation. SPS models ( |Bruzual fc Charlot||2003[ ) 
predict that, following a continuous burst of star 
formation that lasts for ~2x 10 7 yr, the Q(rJgi) 
index increases its value during ~ 2.6 x 10 7 yr and 



then starts to decline (cf., GG96; Martmez-Garcfa 
et al. 2009a)[^] Q(rJgi) mosaics of our sample 



were transformed to polar coordinates rj Spiral 
arms were then "straightened", by adding to the 
9 coordinate a phase as a function of In R, until the 
arms appear horizontal. Unlike our previous stud- 
ies, where the "straightening slope" was different 
depending on the specific arm and region under 
analysis, here we adopt the same "straightening 
slope" for both (or all the) arms in the same ob- 
ject. For each galaxy, this slope is obtained from 
the median of the pitch angles measured in diffcr- 

and Table 



■5.1 



ent wavelengths (see Section 

We divide the "straightened" images in verti- 
cal sections, equivalent to concentric annuli in a 
deprojected image of the galaxy, in order to per- 
form a Fourier analysis^] Here we aim to compare 
the observed Q(rJgi) index profiles with the SPS 
models adopting Fourier techniques. SPS models 
give Q(rJgi) as a function of age, t agc , while ob- 
servations provide Q(rJgi) versus azimuthal an- 
gle (or azimuthal distance after assuming a dis- 
tance to the galaxy). The age gradient signal can 
be thought of as an asymmetric "pulse" (see Fig- 



11 These models were computed for a mixture of young (~ 2% 
by mass) and old stars (with an age of 5 X 10 9 yr). 
Both populations have solar metallicity and a Salpeter ini- 
tial mass function (IMF), with Mi owcr = 0.1 Mq, and 
A^upper = 10 Mq. The choice of M uppcr is justified by the 
observed inverse correlation between the locations where 
color gradients are detected, and those with important H a 
emission (GG96;|Martmez-G arcia et al.||2009a| |Marti'nez-| 



Garcia fe Gonzalez- L6pezrira|2011^ 
12 9 vs. In R; logarithmic arms appear as straight lines in such 
a map, with slope m = cot (— i), where i is the arm pitch 
angle. 

13 We define the median pitch angle for each object by ar- 
ranging the angles in g, r, i, and J from the lowest to the 
highest value, and averaging together the two middle val- 
ues. 

14 We ran some tests that demonstrate that the results are 
independent of the width of the annuli. Widths were chosen 
to be the same as in |Martmez-Garci'a et al.| l |2009a[ | . 



urej5j), whose width is a function of radius; the 
width approaches zero at CR and has larger values 
at smaller radii. The width of this pulse is analo- 
gous to a period, Tq, of a sinusoidal signal. Hence, 
the pulse can also be characterized by its angular 
frequency loq = 2tt/Tq, and analyzed with Fourier 
techniques. 

We compute the Fourier transform of the 
Q(rJgi) index profiles in the annuli, 



f(u> Q )= / Q{e)e-™«d6, 



(2) 



where 



Q(0) 



Q(9) if < 9 < 2tt, 

Q if 9 > 2tt. 



In general, /(loq) is complex, say, of the form 
He(tjQ) + ilm(ojQ). The Fourier spectrum is ob- 
tained as |/(u;q)| = \/Re 2 + Im 2 . We search the 
spectrum for frequencies with amplitudes between 
0.02 and 0.06 mag. This range of amplitudes corre- 
sponds to the ones achieved by the Q(rJgi) index 
of our models)^] including the photometric error 
of the observed Q(rJgi) profiles, ctq ~ 0.05. If 
more than one frequency is found, we then select 
the one with the highest frequency value P'] 

The expected radial behavior of ojq after ap- 
plying the procedure to the MHD simulation 



of Martmez-Garcfa et al. ( 2009b I is shown in Fig- 
ure |6j with a dotted line. We remind the reader 
that the simulation has a spiral pattern with a 
constant angular speed that does not depend on 
radius. The plot also shows (solid line) the behav- 
ior of loq when newborn stars move in perfectly 
circular orbits. In both cases, since the width 
of the pulse goes to zero at the CR radius, its 
reciprocal, uiq, diverges there; Rqr/ Rqlr ~ 0.6. 



15 We consider two models in this analysis (see Figure [sj. 
Model "A" is obtained with the SPS models only. Model 
"B" includes both the dissolution of stellar groups | |Wielen| 

|1977[l afte r 50 Myr, and stellar orbit diffusion. According 
to Widen ( 1977), the diffusion of stellar orbits can enhance 
the dissolution of young stellar groups by increasing their 
internal velocity dispersion. 

16 In general, the lower frequency represents the overall shape 
of the Qlrjgi) profile, whereas higher frequencies corre- 
spond to small-scale structure. 
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Figures [8] through 14 left panels, show the re- 
sults of this analysis when applied to our sample 
of objects. Data uiq values are graphed versus 
-Rmean/-Rend> where i? moa n is the mean galacto- 
centric radius, and -R^n™ ^ s * ne spiral maximum 
radial extent. The latter is determined visually in 
the NIR band indicated in Table [H 

4.2. Azimuthal Phases of the g, r, i, and J 
Bands 



If DWs indeed produce shocks and induce star 
formation, we expect increasingly older stellar 
populations (with intensity peaks at longer wave- 
lengths) to be sequentially located downstream 
(in the direction of rotation) from the shock posi- 
tion within CR, and to be sequentially placed up- 
stream (opposite rotation) outside CR. No gradi- 
ents should be observed at CR. We use also Fourier 
techniques to analyze our sample in a systematic 
manner, at all galactic radii, in search of this ef- 
fect. 

We adopt the method of |Puerari fc DottorT| 
(19971) to determine the phases (basically, the con- 
jugate of the azimuthal positions) of stellar pop- 
ulations of different ages, as traced in the g, r, i, 
and J bands. The method is based on computing 
the Fourier transform of the form 



/(to) 



I R (6)e- im6 d9, 



(3) 



where In is the intensity of radiation, and with 
phase 



$ = tan 



-i J gg[/M] 
\lm[/(m)] 



(4) 



Given that we have selected objects with a strong 
two-armed pattern, m = 2. This means that 2ir 
radians include two periods of the analyzed signal, 
i.e., there is a it radian symmetry. Although this 
may not be exactly true for real arms in spirals 
(P. Grosb0l 2012, private communication), it is a 
very good approximation for our purposes. We an- 
alyze the same concentric annuli in the "straight- 
ened" images as before. 

For the data analysis we assume all spirals in 
our sample trail, and hence that the spiral opens in 
the direction opposite rotation. Furthermore, we 



adopt the convention that the angular coordinate 
increases with rotation, regardless of whether the 
arms have an "S" or a "Z" on-the-sky view. Under 
these premises, for a wave of the form cos((9 + $), 
where <E> is the phase, $ should decrease with wave- 
length inside CR and increase outside CR. There 
should be no phase difference with wavelength at 
CR (see Figure [7]). The results of this exercise are 
shown in Figures [8] through [14} right panels. 

5. Discussion 

For seven objects (NGC 918, NGC 1417, 
NGC 1421, NGC 3001, NGC 3338, NGC 4254, 
and NGC 7125, see Table [1}, we find similari- 
ties between the theoretical expectations and the 
observations, regarding both the radial distribu- 
tion of ujq and the radial run of intensity phases 
with wavelength, at least for some range of radii. 
However, all the objects also show discrepancies. 
For example, in NGC 918, NGC 1417, NGC 1421, 
NGC 3338, NGC 4254, and NGC 7125, the phases 
of all four passbands overlap between w 0.2 and 
0.4 i? mean /i?^™, but then behave in accordance 
to theory up to what is likely the CR radius, as 
judged by the position where luq diverges and 
where, again, phases are the same for all wave- 
lengths. In NGC 918 and NGC 4254, phase dif- 
ferences agree with theory within CR, but then 
do not change direction. This could be due to a 
weaker spiral potential and, hence, weaker shocks 
at radii beyond CR. Conversely, in NGC 3001 
the phases show the expected behavior only out- 
side the CR position, if it is indeed located where 
the luq plot shows a clear divergence. The inten- 
sity phases in NGC 1703 (which may be an SB 
type galaxy) behave as predicted by theory, but 
ujq does not diverge at the CR position expected 
from the phases. 

From this analysis, there is not even one ob- 
ject for which both diagnostics agree perfectly 
with theory at all radii. These findings confirm 
more objectively a previous result from visual in- 
spection: age/color textbook gradients that run 
coherently along entire arms and change direc- 
tion at CR are very hard to detect in practice. 
Possible causes may be a low S/N of the data, 
or masking of the gradients by the intense and 
clumpy (continuum and line) radiation from mas- 
sive star clusters. In the case of the intensity 
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Table 3 

Spiral Radial Maximum Extent 



Galaxy 


-^cnd (arcscc) 




(kpc) 


NGC 


578 


95.7 ± 2.9 (K„) 


10.5 


± 


0.9 


NGC 


918 


75.4 ± 5.8 (-ft"') 


7.9 


± 


0.7 


NGC 


1417 


60.0 ± 5.0 (J) 


16.6 


± 


1.4 


NGC 


1421 


92.8 ± 2.9 (K e ) 


13.2 


± 


1.1 


NGC 


1703 


27.5 ± 1.4 (H) 


2.7 


± 


0.3 


NGC 


3001 


44.0 ± 2.8 (H) 


7.6 


± 


0.8 


NGC 


3338 


145.0 ± 5.0 (K s ) 


16.7 


± 


1.6 


NGC 


4254 


157.5 ± 7.5 (K s ) 


12.6 


± 


0.8 


NGC 


4603 


45.4 ± 1.4 (H) 


6.5 


± 


0.6 


NGC 


4939 


145.0 ± 5.0 (K s ) 


32.7 


± 


2.8 


NGC 


6907 


55.1 ± 2.9 (K s ) 


12.7 


± 


1.3 


NGC 


6951 


105.0 ± 5.0 (K s ) 


12.7 


± 


1.1 


NGC 


7125 


95.7 ± 2.9 (J) 


20.7 


± 


1.8 



Note. — Columns 2 and 3: spiral radial maxi- 
mum extent (determined visually), in arcscc and 
kpc. respectively, and bandpass used to determine 
it. 



phases method, dusty environments can attenuate 
the shorter wavelengths, hindering the detection of 
an ordered age sequence^ It is also possible that 
the gradients are perturbed by the "infant mortal- 
ity" of star clusters ( Lada fc Lada|2003 ), or by the 
destruction of an ordered age front by supernovae 
shocks and stellar winds ( Mueller fc Arnett||1976 



the IGIM F (the Integrated Galaxial IGIMF) 
theory of Pflamm-Altenburg & Kroupa (20081 



Gerola & Seiden 1978). On the other hand, even 



in the presence of DWs, it is possible that no gra- 
dients arc produced. For example, there may be 
inadequate physical conditions for the creation of 
the gradients, related to the strengths of the shock 
or of the magnetic field ( Efremov||2010 ). There is, 
of course, the possibility that a real constant pat- 
tern speed for all radii does not exist, i.e., that 
the spiral structure is not explained by the DW 
theory. Our subsample, however, was selected ex- 
plicitly to minimize this eventuality. Furthermore, 
we had previously established the link of the gra- 
dients in the sample, taken as an ensemble, to disk 
dynamics. As confirmed by the comparison with 
an MHD simulation with a constant spiral pat- 
tern speed, an apparent radial dependency of Q p 
is caused by neglecting non-circular streaming mo- 
tions in the analysis of the gradients. 

Conversely, there could be color gradients 
that are not linked to the dynamics of the 
disk. One speculative scenario is provided by 



and Pflamm-Altenburg et al. (2009). This theory 



is able to explain the Ha cutoff in galaxy disks 
without a corresponding cutoff in star formation 
(as deduced from the non-ionizing far-ultraviolet 
flux given out by B-type stars) . The IGIMF theory 
posits that maximum star cluster mass depends on 
gas surface density (or, equivalently, on star for- 
mation rate){^] and that maximum stellar mass 
depends on cluster mass. If this is indeed the case, 
one could obtain a color gradient from coeval star 
clusters across a region with a gas surface density 
gradient. 

At any rate, a statistical analysis, tied to spiral 
dynamics, of color gradients in several objects is 
necessary to identify the gradients that are linked 
to DWs. This is even more important since by 
general rule gradients are painstakingly identified 
only in selected regions and not along entire arms. 

It is, however, important to stress that our new 
results imply that long-lived modes seem to have 
produced gradients that run for a sizeable part of a 
spiral arm in at least 50% of the analyzed objects. 
This is, in itself, remarkable. 

For the rest of the objects (NGC 578, NGC 1703, 
NGC 4603, NGC 4939, NGC 6907, and NGC 6951), 



17 This was the original motivation behind the development 
of the reddening-insensitive index Q(rjgi). 



18 This has been corroborated by Gonzalez-Lopezlira et al. 
(2012) in M33. 
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the Fourier methods failed to detect signs of gradi- 
ents besides the ones that had already been identi- 
fied by eye. It is interesting to point out, though, 
that NGC 4603 shows a multi-armed structure in 
the optical and NIR bands, although it was kept 
in the sample because of its two-armed structure 
in the mass map. For NGC 6951 (type SAB in the 
RC3) and NGC 6907 (type SB in the RC3), the 
CR radius may lie around -Rgjj? ~ 0.8. This is 
possibly an indication of color gradients within the 
bar region, as already found in Martrnez-Garcfa 



& Gonzalez-Lopezlira (2011) 



We find no relation between Hubble type and 
success rate in our Fourier tests. 

5.1. Pitch Angles Test 



As a complementary test of DW presence, we 
examine the pitch angles of the spiral arms in dif- 
ferent wavebands (g, r, i, and J), and in the "re- 
solved mass maps" previously obtained with the 
method of Zibetti et al. ( 2009 1 . As already men- 



tioned in the introduction, according to Efrcmov 
(2010), age gradients are observed across arm seg- 



ments that fulfill the conditions for a strong shock 
to occur. These segments seem to have an irreg- 
ular magnetic field, while the regions with a reg- 
ular magnetic field harbor instead weak shocks, 
and chains of star complexes with a nearly regu- 
lar spacing. Strong shocks are associated with a 
larger pitch angle of the spiral segment. 

Another prediction involving pitch angles in the 
DW theory paradigm comes from the local disper- 
sion relation. At a fixed mass surface density, the 
pitch angle increases proportionally to the square 



of the velocity dispersion (see, e.g., 


Hozumi||2003 


Athanassoula et al. 


2010 Martmez-Garcfa 


2012), 



that is, with the age of the stellar population and, 
hence, with image wavelength in the optical and 
NIR. 

Pitch angles were measured with Fourier tech- 



niques (Considere & Athanassoula 1982| 1988 


Puerari &; Dottori 


1992 


Saraiva Schroeder et al.| 


1994 


Davis et al. 


2012 


Savchenko 


2012 


), in the 



radial range where the spiral perturbation seems 
more prominent (see Table [5| . This method as- 
sumes a logarithmic geometry for the spirals, al- 
though this might not be the case for all objects 
at all radii ( Ringermacher & Mead 2009 ) . 



Figure 15 shows histograms of P x - Pj, where P 
is the pitch angle and x can be g, r, i, or "mass", 
depending on the image where P has been mea- 
sured^ Negative values of the differences indicate 
that the pitch angle is larger in the NIR than in 
the optical. In the histograms of P g — Pj, P r — Pj, 
and Pi — Pj, there are indeed more objects with 
negative differences. The median values of each 
histogram are, respectively, -0.9°, -1.1°, and 0.0°. 
This result is in agreement with |Grosb0l &: Pat-| 



SIS 



( 1998 ), who find tighter arms in bluer colors in 
images of four spirals, suggesting the presence of 
DWs. 

Interestingly, the median value of the P maS s — 
Pj histogram is 2.2°, so that the mass arms are 
statistically more open than the arms in NIR light. 
An inspection of the mass maps indicates that the 
spiral arms loci in the NIR images do not always 
coincide with the location of the arms in the mass 
maps (see Figures [3] and a fact that can be 
explained by the contribution of young stars to 
the NIR light. 

Finally, we compare the median pitch angle dis- 
tribution in those galaxies that we have classified 
as having radially extended azimuthal color gradi- 
ents, against the distribution in the galaxies where 
we have detected localized gradients only. Fig- 
ure [16] displays the pitch angle histograms for each 
data set. The statistical significance is low, due to 
the small number of objects, but the objects with 
signs of widespread color gradients may have more 
open spiral arms (larger pitch angles). 

6. Conclusions 



From a Fourier analysis of 13 objects, with a 
spiral arm structure that can be likely explained 
by DW theory, we find that ~ 50% of them show 
evidence of color gradients that run for at least a 
range of radii. For the remaining objects, the ev- 
idence of extended gradients is insufficient. This 



notwithstanding, our previous studies (Gonzalez 



lf Absolute values from Table 5 were adopted for analysis 



However, P x values were discarded if their sign does not 



agree with that of Pj. (Signs indicate whether the galaxy 
has an "S" or a "Z" on-the-sky view.) The pitch angles 



measured for NGC 578 with the Fourier method were also 



|rejected, since their nominal sign is at odds with the ob- 
served on-the-sky view. 
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Table 4 



Corotation Positions (-Rcr) 





Galaxy 




Ren (arcsec) 


Rcr (kpc) 


NGC 918 


0.63 


47.5 ± 3.7 


5.0 ± 0.4 


NGC 1417 


0.67 


40.2 ± 3.4 


11.1 ± 0.9 


NGC 1421 


0.49 


45.5 ± 1.4 


6.5 ± 0.5 


NGC 3001 


0.55 


24.2 ± 1.5 


4.2 ± 0.4 


NGC 3338 


0.38 


55.1 ± 1.9 


6.3 ± 0.6 


NGC 4254 


0.46 


72.5 ± 3.5 


5.8 ± 0.4 


NGC 7125 


0.67 


64.1 ± 1.9 


13.9 ± 1.2 



Note. — Column 1: object with signatures of a radially ex- 
tended azimuthal gradient. Column 2: corotation likely po- 
sition, fl mC an/-Rcnd ■ The mean value for the six objects is 
flmean/fiend — 0.55. Columns 3, and 4: corotation radius, 
RcH: m arcscconds and kiloparsccs, respectively. 



I Table 5 1 

I Pitch Angle Values ~l 



Galaxy 


A/2(arcsec) 


P, 






Pr 




P 






Pj 




Pm 






"Median" 


NGC 578 


(20.3-63.8) 


-34.70 


m- 


-2) 


-35.72 


[ra 


=2) 


-37.92 


ra 


=2) 


-41.68 


[ra 


=2) 


-49.45 


[ra 


=2) 


40.97 a 


NGC 918 


(18.6-57.2) 


20.60 


m- 


2) 


21.77 


[ra 


=2) 


22.63 


ra- 


= 2) 


23.54 


[ra 


=2) 


29.42 


[ra 


=2) 


22.20 


NGC 1417 


(20.0-60.0) 


36.82 


m= 


=2) 


33.76 


m 


=2) 


36.82 


rn 


=2) 


41.71 


[ra 


=2) 


39.13 


[in 


=2) 


36.82 


NGC 1421 


(34.4-91.6) 


36.17 


m- 


2) 


37.17 


[ra 


=2) 


38.21 


va- 


= 2) 


39.31 


[ra 


=2) 


-81.66 


[ra 


=2) 


37.69 


NGC 1703 


(12.4-27.5) 


16.55 


m= 


,2) 


16.29 


m 


=2) 


15.79 


ra 


=2) 


-36.16 


[ra 


=1) 


-55.62 


[in 


=2) 


16.04 


NGC 3001 


(19.2-44.0) 


31.55 


m= 


2) 


31.55 


m 


=2) 


32.45 


ra 


=2) 


32.45 


[ra 


=2) 


34.40 


[m 


=2) 


32.00 


NGC 3338 


(45.0-145.0) 


-13.43 


m- 


=2) 


-13.60 


[ra 


=2) 


-13.77 


ra 


= 2) 


-13.77 


[ra 


=2) 


83.87 


[ra 


=2) 


-13.68 


NGC 4254 


(30.0-157.5) 


-29.30 


m= 


=3) 


18.13 


m 


= 1) 


-29.30 


ra 


=3) 


-29.30 


[ra 


=3) 


22.62 


[m 


=1) 


-29.30 


NGC 4603 


(13.8-45.4) 


25.61 


m- 


2) 


26.87 


[ra 


=2) 


27.54 


ra 


= 2) 


28.25 


[ra 


=2) 


30.59 


[ra 


=2) 


27.20 


NGC 4939 


(50.0-145.0) 


10.96 


m= 


,2) 


11.07 


m 


=2) 


11.07 


ra 


=2) 


10.74 


[ra 


=2) 


29.89 


[in 


=3) 


11.02 


NGC 6907 


(34.8-55.1) 


24.89 


m= 


:2) 


24.34 


m 


=2) 


27.34 


ra 


=2) 


27.34 


[ra 


=2) 


45.15 


[in 


=2) 


26.12 


NGC 6951 


(45.0-105.0) 


-31.48 


m- 


2) 


-33.26 


[ra 


=2) 


-32.35 


ra- 


= 2) 


-34.23 


[ra 


=2) 


-36.30 


[ra 


=2) 


-32.80 


NGC 7125 


(20.3-95.7) 


-23.88 


m= 


=2) 


-26.76 


m 


=2) 


-35.99 


ra 


=2) 


-40.84 


[ra 


=2) 


-52.36 


[in 


=2) 


-31.38 



NOTE. — Column 1: analyzed object. Column 2: radial ranges where the spiral seems more prominent, in arcscconds. Columns 
3-6: pitch angles, in degrees, for the g, r, i, and J bands, respectively. Negative values indicate a "Z" on-the-sky view for the 
object, while positive values indicate an "S" on-thc-sky view. Absolute values were adopted for analysis. The mean error is ~ 1° 
(without taking into account the error due to the deprojection parameters). The dominant Fourier mode, m, is indi cated in 
parenthe sis. Column 7: pitch angle, in degrees, for the resolved maps of stellar mass surface density ("mass map" , [Zibetti etj 
|al.|2009| . Column 8: median pitch angle, in degrees, adopted to "straighten" the spiral arms (see Section |4.1| l. We define the 
median pitch angle for each object by arranging the P g , P r , Pi, and Pj values from lowest to highest, and averaging together 
the two middle values. 

a Computed by fitting the arms' slope in the 8 vs. In R polar image (sec, e.g., | Martinez- Garcia | |2012| l. The Four icr method 
yields a "Z" on-thc-sky view, but the object actually has an "S" view. 
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fe Graham||l996l |Martmez-Garcia et al.||2009a|b 



Martmez-Garcia fe Gorizalez-LopezlirapOlT ) sug- 
gest that azimuthal age/color gradients (or candi- 
dates) can be found for most galaxies, at least in 
some regions of the spiral arms. A relation of such 
very localized gradients with the dynamics of the 
disk has been previously established by the com- 
parison with an MHD simulation with a constant 



pattern speed at all radii (Martmez-Garcia et al 



2009b) 



Our results are consistent with other observa- 
tional studies of age patterns (or offsets) across 



spiral arms 1 


Efremov 


1985 Tamburro et al. 


2008 


Egusa et al. 


2009 


Grosb0l & Dottori 2009| 


Sanchez-Gil et al. 


2011) 


. By means of an Ha to 



far-UV flux ratio method, Sanchez-Gil et al. (2011 ) 



found age gradients across the spiral arms of the 
grand-design spirals M74 and M100 (NGC 628 and 
NGC 4321, respectively). Likewise, the "breaks" 
in the radial metallicity distribution found near 
CR in a sample of ~ 20 galaxies by Scarano "et| 



al. (2011) and Scarano & Lepine (20131 imply 
that spiral arms must be long-lived structures. A 
dominant pattern speed must exist with a unique 
CR radius, or otherwise any discontinuities in 
the radial metallicity profiles would be smoothed 
out ( [Scarano fe Lepine|[2013| . 



On the other hand, simulations have shown mu- 
tually exclusive results. Long-lived quasi-steady 
features have been obtained for some models (e.g., 



Donner & Thomasson 


1994 


Zhang 


1998 


I. Other 


works, however (e.g., Sellwood||2011 


Wada et al. 


2011 


Fujii et al.|2011 ) 


, have found a lack of a per- 



sistcnt pattern, and instead recurrent short-lived 
transient spirals only. 

Any model that intends to explain spiral struc- 
ture has to avoid the winding dilemma, that ul- 
timately led to the idea of stationary DWs with 
a fixed pattern speed for all radii. For this same 
reason, the observational results of |Merrifield et| 



al. (20061 and Meidt et al. (2009) for some ob- 



jects, where fl p decreases with radius, imply the 
existence of transient spiral modes that emerge 
and disappear at different stages in the disk evolu- 
tion. Such transient arms cannot produce stellar 
age/color gradients, though. 

A theory that reconciles the existence of 
age/color gradients with a radially varying pat- 
tern speed is still missing. However, as argued 
in Section [T]T] the observation of a pattern speed 



that varies with radius may be an artifact pro- 
duced by the non-circular streaming motions of 
the stars newly born in shocked material. This 
effect can explain the coexistence of azimuthal 
age/color gradients across spiral arms, and of a 
spiral pattern speed f2 p that seems to vary with 
radius in the same way as the orbital frequency tt. 

We acknowledge Dr. Yuri Efremov, the ref- 
eree, for his comments and suggestions that have 
greatly improved this paper. Special thanks go 
to Ivanio Puerari, for helpful discussions about 
Fourier analysis of spiral galaxies. E.M.G. ac- 
knowledges a CONACYT grant for a postdoctoral 
fellowship at INAOE, and former postdoctoral fi- 
nancial support from UNAM (DGAPA), Mexico. 
R.A.G.L. is also grateful for financial support from 
both CONACYT and D G APA-UN AM . 

REFERENCES 

Arzoumanian, D., Andre, P., Didelon, P., et al. 

2011, A&A, 529, L6 

Athanassoula, E., Romero-Gomez, M., Bosma, A., 
& Masdemont, J. J. 2010, MNRAS, 407, 1433 

Block, D. L., Bertin, G., Stockton, A., et al. 1994, 
A&A, 288, 365 

Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 
1000 

Clarke, O, & Gittins, D. 2006, MNRAS, 371, 530 

Chariot, S., & Fall, S. M. 2000, ApJ, 539, 718 

Considere, S., & Athanassoula, E. 1982, A&A, 
111, 28 

Considere, S., & Athanassoula, E. 1988, A&AS, 
76, 365 

da Cunha, E., Chariot, S., & Elbaz, D. 2008, MN- 
RAS, 388, 1595 

Davis, B. L., Berrier, J. O, Shields, D. W., et al. 

2012, ApJS, 199, 33 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, 
H. G., Jr., et al. 1991, Third Reference Cat- 
alogue of Bright Galaxies. (Berlin: Springer), 
(RC3) 



11 



Dobbs, C. L., & Bonnell, I. A. 2008, MNRAS, 385, 
1893 

Dobbs, C. L., & Price, D. J. 2008, MNRAS, 383, 
497 

Donncr, K. J., & Thomasson, M. 1994, A&A, 290, 
785 

Efremov, Y. N. 1985, Soviet Astronomy Letters, 
11, 69 

Efremov, Y. N. 2010, MNRAS, 405, 1531 

Efremov, Y. N. 2011, Astronomy Reports, 55, 108 

Egusa, F., Kohno, K., Sofue, Y., Nakanishi, H., & 
Komugi, S. 2009, ApJ, 697, 1870 

Ferreras, I., Cropper, M., Kawata, D., Page, M., 
& Hoversten, E. A. 2012, MNRAS, 424, 1636 

Foylc, K., Rix, H.-W., Dobbs, C. L., Leroy, A. K., 
& Walter, F. 2011, ApJ, 735, 101 

Fujii, M. S., Baba, J., Saitoh, T. R., et al. 2011, 
ApJ, 730, 109 

Gerola, H., & Seiden, P. E. 1978, ApJ, 223, 129 

Gittins, D. M., & Clarke, C. J. 2004, MNRAS, 
349, 909 

Goodwin, S. P., & Bastian, N. 2006, MNRAS, 373, 
752 

Gonzalez, R. A., & Graham, J. R. 1996, ApJ, 460, 
651 (GG96) 

Gonzalez-Lopezlira, R. A., Pflamm-Altenburg, J., 
& Kroupa, P. 2012, ApJ, 761, 124 

Grand, R. J. J., Kawata, D., & Cropper, M. 2012a, 
MNRAS, 421, 1529 

Grand, R. J. J., Kawata, D., & Cropper, M. 2012b, 
MNRAS, 426, 167 

Grosb0l, P. J., & Patsis, P. A. 1998, A&A, 336, 
840 

Grosb0l, P., & Dottori, H. 2009, A&A, 499, L21 

Grosb0l, P., Dottori, H., & Grcdel, R. 2006, A&A, 
453, L25 



Hozumi, S. 2003, in Galaxies and Chaos, ed. G. 
Contopoulos & N. Voglis (Lecture Notes in 
Physics, Vol. 626; Berlin: Springer), 380 

James, P. A., & Seigar, M. S. 1999, A&A, 350, 791 

Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 

Lin, C. C, & Shu, F. H. 1964, ApJ, 140, 646 

Lindblad, B. 1963, Stockholms Observatoriums 
Annaler, 22, 5 

Martmez-Garci'a, E. E., Gonzalez-Lopezlira, 
R. A., & Bruzual-A, G. 2009a, ApJ, 694, 512 

Martmez-Garci'a, E. E., Gonzalez-Lopezlira, 
R. A., & Gomez, G. C. 2009b, ApJ, 707, 1650 

Martmez-Garci'a, E. E., & Gonzalez-Lopezlira, 
R. A. 2011, ApJ, 734, 122 

Martmez-Garci'a, E. E. 2012, ApJ, 744, 92 

Meidt, S. E., Rand, R. J., & Merrifield, M. R. 
2009, ApJ, 702, 277 

Merrifield, M. R., Rand, R. J., & Meidt, S. E. 
2006, MNRAS, 366, L17 

Mould, J. R. ct al. 2000, ApJ, 529, 786 

Mueller, M. W., & Arnett, W. D. 1976, ApJ, 210, 
670 

Patsis, P. A., Heraudcau, P., & Grosb0l, P. 2001, 
A&A, 370, 875 

Paturel, G., Fang, Y., Petit, C, Garnier, R., & 
Rousseau, J. 2000, A&AS, 146, 19 

Paturel, G., Theureau, G., Bottinelli, L., 
Gougucnhcim, L., Coudreau-Durand, N., Hal- 
let, N., Petit, C. 2003, A&A, 412, 57 

Pflamm-Altenburg, J., & Kroupa, P. 2008, Na- 
ture, 455, 641 

Pflamm-Altenburg, J., Weidner, C, & Kroupa, P. 
2009, MNRAS, 395, 394 

Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et 
al. 2010, A&A, 518, LI 

Puerari, I., & Dottori, H. A. 1992, A&AS, 93, 469 

Puerari, I., & Dottori, H. 1997, ApJ, 476, L73 



12 



Rix, H.W., & Rieke, M. J. 1993, ApJ, 418, 123 

Rhoads, J. E. 1998, AJ, 115, 472 

Roberts, W. W. 1969, ApJ, 158, 123 

Roca-Fabrega, S., Valenzuela, O., Figueras, F., et 
al. 2013, MNRAS, submitted 

Ringermachcr, H. I., & Mead, L. R. 2009, MN- 
RAS, 397, 164 

Sanchez-Gil, M. C, Jones, D. H., Perez, E., et al. 
2011, MNRAS, 415, 753 

Saraiva Schroeder, M. F., Pastoriza, M. G., Ke- 
pler, S. O., & Puerari, I. 1994, A&AS, 108, 41 

Savchenko, S. S. 2012, Astrophysical Bulletin, 67, 
310 

Scarano, S., Jr., Lepine, J. R. D., & Marcon- 
Uchida, M. M. 2011, MNRAS, 412, 1741 

Scarano, S., & Lepine, J. R. D. 2013, MNRAS, 
428, 625 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 
1998, ApJ, 500, 525 

Sellwood, J. A. 2011, MNRAS, 410, 1637 

Speights, J. C, & Westpfahl, D. J. 2012, ApJ, 752, 
52 

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 
2006, AJ, 131, 1163 

Tamburro, D., Rix, H.-W., Walter, F., et al. 2008, 
AJ, 136, 2872 

Thuan, T. X., Gunn, J. E. 1976, PASP, 88, 543 

Wada, K., Baba, J., & Saitoh, T. R. 2011, ApJ, 
735, 1 

Wade, R. A., Hoessel, J. G., Elias, J. H., Huchra, 
J. P. 1979, PASP, 91, 35 

Wielen, R. 1977, A&A, 60, 263 

Zhang, X. 1998, ApJ, 499, 93 

Zibetti, S., Chariot, S., & Rix, H.-W. 2009, MN- 
RAS, 400, 1181 

Zibetti, S. 2009, |arXiv:0911.4956| 



This 2-column preprint was prepared with the AAS IATj^X 
macros v5.2. 



13 



50 - 



40 - 



"i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



Input Q 
Output Q p 

Orbital frequency Q 




6 



7 8 
R (kpc) 



9 



10 



Fig. 1. — Pattern speed, 17 p (km s -1 kpc -1 ), vs. radius, R (kpc), obtained from the MHD simulation 
analyzed in Martmez-Garci'a et al. ( |2009b ). Long-dashed line: input Sl p ; solid line: pattern speed values 
obtained by applying the GG96 method to synthetic azimuthal color(age) gradients across the arms, under 
the circular motion assumption; dotted line: orbital frequency Q. 
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Fig. 2. — Mass-to-light ratios (M/Ls) as a lunction of the color-color space (g — i) vs. (i — J). Left panel: 
median logarithmic M/L in the J band, (M/L)j, for models binned in elements with size 0.05 x 0.05 mag 2 
in two-dimensional color. Right panel: rms of log M/L. M/L is the ratio between stellar mass and the 
light that reaches the observer, i.e., the effective M/L (Zibetti et al. 2009). Optical g and i photometry 
in the Thuan-Gunn system ([Thuan fc Gunn||1976{ |Wade et al.||1979[ ). J band calibrated as in the 2MASS 
survey ( |Skrutskie et~aL|2006[ ) . 
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Fig. 8. — Panels (a) and (c): frequencies ojq vs. ^mcan/^™ f° r NGC 578, and NGC 918, respectively. 
Panels (b), and (d): two-armed phases (see Section 4.2), in the g (green solid line), r (red dotted line), i (blue 
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Fig. 9. — Panels (a), and (c): frequencies ujq vs. radius for NGC 1417 and NGC 1421, respectively. Panels 
(b) and (d): two-armed phases vs. radius for NGC 1417 and NGC 1421, respectively. Corotation likely 
radial positions are indicated. Symbols as in Figure [HJ 
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Fig. 10. — Panels (a) and (c): frequencies ujq vs. radius for NGC 1703 and NGC 3001, respectively. Panels 
(b) and (d): two-armed phases vs. radius for NGC 1703 and NGC 3001, respectively. Corotation likely 
radial position is indicated for NGC 3001. Symbols as in Figure [HJ 
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Fig. 11. — Panels (a) and (c): frequencies ujq vs. radius for NGC 3338 and NGC 4254, respectively. Panels 
(b) and (d): two-armed phases vs. radius for NGC 3338 and NGC 4254, respectively. Corotation likely 
radial positions are indicated. Symbols as in Figure [HJ 
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Fig. 12. — Panels (a) and (c): frequencies ujq vs. radius for NGC 4603 and NGC 4939, respectively. Panels 
(b) and (d): two-armed phases vs. radius for NGC 4603 and NGC 4939, respectively. Symbols as in Figure|8] 
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Fig. 13. — Panels (a) and (c): frequencies ujq vs. radius for NGC 6907 and NGC 6951, respectively. Panels 
(b) and (d): two-armed phases vs. radius for NGC 6907 and NGC 6951, respectively. Symbols as in Figure|8] 
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Fig. 14. — Panel (a): frequencies ujq vs. radius for NGC 7125. Panel (b): two-armed phases vs. radius for 
NGC 7125. Corotation likely radial position is indicated. Symbols as in Figure [HJ 
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Fig. 15. — Histograms of the pitch angle differences between the optical bands g (top left), r (top right), i 
(bottom left), the "mass map" (bottom right), and the NIR J band. 
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Fig. 16. — Histograms of the median pitch angle values. Top: objects with signs of widespread color 
gradients; bottom: objects with localized color gradients only. 
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